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We have discovered 21 millisecond pulsars (MSPs) in the globular cluster Terzan 5 
using the Green Bank Telescope, bringing the total of known MSPs in Terzan 5 
to 24. These discoveries confirm fundamental predictions of globular cluster 
and binary system evolution. Thirteen of the new MSPs are in binaries, of 
which two show eclipses and two have highly eccentric orbits. The relativis- 
tic periastron advance for the two eccentric systems indicates that at least one 
of these pulsars has a mass >1.68Mo at 95% confidence. Such large neutron 
star masses constrain the equation of state of matter at or beyond the nuclear 
equilibrium density. 

The extremely high stellar densities (10^— lO^pc"'^) in the cores of globular clusters (GCs) 
result in stellar interactions which produce and destroy binary systems as well as exchange their 
members ©. Formed by the death of massive stars early in a cluster's history, neutron stars 
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(NSs) usually reside near the cores of clusters due to their relatively large masses and the mass 
segregation induced by dynamical friction. There they are likely to interact with one or more 
stars over the 10^° yr lifetimes of GCs. These interactions lead to a production rate of low-mass 
x-ray binaries (LMXBs) and their progeny such as MSPs [via the recycling mechanism (T)], that 
is highly enhanced compared to the rate in the Milky Way Galaxy. Before our observations, there 
were 80 pulsars (most of them binary MSPs) known in 24 GCs (3), with the relatively massive 
and nearby cluster 47 Tucanae containing 22 of these (Jj. Finding and monitoring many pulsars 
in a single cluster provides unique probes into a range of GC, binary evolutionary, and stellar 
astrophysics In addition, GCs produce exotic systems such as highly eccentric (0) and 

MSP— main sequence star binaries ® . 

The interaction rate between NSs and other stars or binaries in a GC is a complex function of 
the total cluster mass, the size of the cluster core (O and its stellar density, the initial stellar mass 
function, and the level of mass segregation present in the cluster tiiQi) . However, relatively simple 
theoretical modelling of stellar interaction rates (11), as well as the known LMXB and several 
additional x-ray sources detected recently with the Chandra X-ray Observatory { 12 ) indicate that 
the dense, massive, and metal-rich GC Terzan 5 has one of the highest stellar interaction rates 
of any cluster in the Galaxy dTib and perhaps also the largest number of MSPs (|7?1I75|) . But 
because Terzan 5 [Galactic coordinates (/,6)=(3?8, 1?7)] is distant (Z}=8.7±2kpc), and located 
within ~1 kpc of the Galactic center (176b . the large column density of interstellar free electrons 
(i.e. the dispersion measure, DM~240pccm~^) produces considerable dispersive smearing (oc 

for typical pulsar search data) and scatter broadening (oc u^^'*) of radio pulses, which hinders 
pulsation searches for MSPs at the often observed radio frequencies of z/=400— 1400 MHz. 

Deep radio images of Terzan 5, which are not affected by the dispersive effects of the inter- 
stellar medium (ISM), were made with the Very Large Array (VLA) and showed what appeared 
to be the integrated emission from many tens or even hundreds of pulsars with low flux densi- 
ties (|77b . Surprisingly, though, numerous deep searches for pulsars over the past 15 years using 
the Parkes radio telescope at 1400 MHz have identified only three (17^ : the 11 -ms Ter 5 A with 
its unusual eclipses dTPb . and the two isolated MSPs Ter 5 C (1201) and D (1^ . 

Observations and Data Analysis On 17 July 2004, we observed Terzan 5 for 5.9 hr using 
the National Radio Astronomy Observatory's (l22i) 100-m Green Bank Telescope (GBT) with 
the Pulsar Spigot backend IQS^. The S-band receiver provided 600-MHz (1650-2250 MHz) of 
relatively interference-free bandwidth in two orthogonal polarizations, which the Spigot summed 
and synthesized into 768x0.78125 MHz frequency channels every 81.92 /xs. With the known 
DM and scattering time scale ( (2?| l toward the cluster, the effective time resolution of the data 
was ~0.3 ms. We observed only a single position since the 6^5 telescope beam width at these 
frequencies is significantly larger than the 0^83 half-mass radius of the cluster (l25b . The large 
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unblocked aperture of the GBT, the wide bandwidth provided by the Spigot and S-band receiver, 
and the move to higher observing frequencies together provided an increase in sensitivity over 
the Parkes searches (B71) by factors of ~5 for typical recycled pulsars and by >10 for MSPs with 
spin periods Ppsr^2 ms (l26t . 

We searched the observation by de-dispersing the raw data into 40 separate time series with 
DMs ranging from 230— 250pccm^^ and spaced by O.Spccm^'^. We Fourier-transformed the 
full 5.9-hr time series as well as 10-, 20- and 60-min sections and searched them using Fourier- 
domain acceleration search techniques STTh in a manner similar to that described in Ransom 
et al. in order to maintain sensitivity to pulsars in compact binary systems. In that first 
observation we discovered 14 new pulsars, Ter 5 E— R (OPb . Using eight more observations taken 
between July and November 2004, we found seven additional pulsars and determined the basic 
orbital parameters of 10 of the 13 new binaries (Fig. 1 and Table 1). In general, the uncalibrated 
flux densities of the pulsars (as computed by comparing the integrated signal from a pulsar to 
the predicted total system noise level) were constant to within a level of <50% between these 
observations implying that diffractive scintillation will not largely affect pulsar measurements 
for Terzan 5 at these observing frequencies. Several additional binary pulsar candidates from 
these data remain unconfirmed but may reappear in future observations at more fortuitous orbital 
phases. Since the volume enclosed by the GBT beam at 1950MHz out to a distance of lOkpc 
is ~10^'^ kpc'^ and it has been estimated that there are roughly 10^ observable MSPs within 
the Galaxy (with volume ~10^ — 10^ kpc^), it is possible that a detectable foreground MSP is 
within our beam. However, given that the DM towards Terzan 5 is known to within ~5%, the 
probability of that MSP having a similar DM is much smaller, and hence we are confident that 
all of the new pulsars are members of Terzan 5. Positions with <1" accuracy of the new MSPs 
from pulsar timing will solidify the associations. 

The New Pulsar Population With our initial discovery data, the ensemble of pulsars all have 
measured DMs with errors <0.1 pccm"'^, flux densities 5*1950 with fractional errors of ~30%, 
and precisely determined spin periods Ppsr. These measurements allow us to compare the new 
systems with known MSPs such as those in 47 Tuc. Once precise positions and pulsar spin period 
derivatives have been determined from timing observations over the next year (the positions of 
the new pulsars are currently known only to within the 6. '5 primary beam of the GBT), a variety of 
tests of the gravitational potential and internal dynamics of Terzan 5 will be possible, in addition 
to estimates of the ages and magnetic field strengths of many of the pulsars (31 ). 

The DMs for the pulsars in Terzan 5 are distributed with a central Gaussian-like component 
with an average and standard deviation of 237.8 ± 1.4pccm~^, and five outliers having higher 
and lower DMs (Table 1). The overall spread in DM, 9.5 pccm"^ (13 times larger than what 
is observed for 47 Tuc), is the largest known for any GC and is likely due to the large average 
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DM for Terzan 5 and irregularities in the ISM along the differing sight- lines toward its pulsars. 
The Gaussian-like component of the DM distribution likely corresponds to a group of centrally 
concentrated pulsars near the cluster core. The position of Ter 5 C, which has an average DM, is 
about 10" or one core radius north of the cluster center and supports such a notion. Likewise, 
pulsars with outlying DMs, such as Ter 5 D and J, are probably more offset like Ter 5 A, which 
has a high DM and is located 36" from the center SWk . 

The VLA radio imaging of Terzan 5(77) revealed numerous point sources within 30" of the 
cluster center, as well as ~2 mJy of diffuse emission at 1400 MHz within ~10". The central emis- 
sion was attributed to 60—200 unresolved pulsars assuming D=7.1kpc, a standard luminosity 
distribution, and a minimum pulsar luminosity (L1400 = 5'i4oo x D'^) of i^i4oo,min=0.3 mJy kpc^. 
A distance to Terzan 5 of D=8.7kpc (I76i) . would increase the number of unresolved pulsars by 
~50%. For a typical pulsar radio spectral index of —1.6 (l32b . and omitting the flux densities 
of Ter 5 A and C which are located outside the region of diffuse emission {Hh, the integrated 
measured flux density from the other 22 pulsars at 1400 MHz is ~1.3mJy. If the five pulsars 
with outlying DMs reside outside the cluster core, the total flux density from the remaining pul- 
sars is only ~1 mJy. None of the new pulsars can individually (or in combination with nearby 
Ter 5 C) account for the bright (1.42 mJy) point source "N" located 12" north of the cluster 
center (ITTb . However, given its very wide pulse profile, lack of a flat off-pulse baseline, and 
positional coincidence (l20t . Ter 5 C could possibly emit significant unpulsed emission like the 
MSP J02 18+4232 (EU and thereby account for aU of the flux density of source N. 

The differential luminosity distribution of the new pulsars resembles that of the 47 Tuc (0), 
M15 (5), and the Galactic disk (34) populations by following the normal dlogN = —dlogL 
relation, with no sign of a downturn at the low luminosity end (Li4oo~l— 2 mJy kpc^). This 
implies that we have not reached the lowest luminosities of the intrinsic pulsar distribution in 
Terzan 5, nor have we yet reached the sensitivity limit of the observing system 05h. It is likely 
that tens of less-luminous pulsars remain to be discovered in Terzan 5, and possibly even several 
bright ones (e.g. the source N). Such bright pulsars may remain undetected if they are members 
of very compact or massive binaries due to eclipsing or excessive Doppler accelerations or if they 
have extremely fast spin periods (Ppsr^l.5 ms). Detailed searches using more advanced search 
techniques may yet uncover these systems. 

The 24 known pulsars in Terzan 5 and the 22 in 47 Tuc appear to have different spin period 
distributions. The 47 Tuc pulsars are a homogeneous population with periods 2.1— 7.6 ms 
while those in Terzan 5 have a flatter distribution that includes six pulsars slower than 7.6 ms 
and the four fastest pulsars known in GCs (Ppsr=1.67, 1.73, 2.05, and 2.07 ms). A Kolmogorov- 
Smimov test suggests that these samples are drawn from different parent distributions at 85% 
confidence. These differences may be related to the dynamical states of the cluster cores [i.e. 
just pre- or just post-core collapse for Terzan 5 (I76b l. different epochs of MSP creation, or the 
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occurrence of unusual evolutionary mechanisms that only manifest themselves at very high (i.e. 
>10^pc^^) stellar densities. 

Individual Pulsars The binaries Ter 5 E and W have orbital periods of Porb~60days and 
~4.9 days, respectively. In the dense stellar environment of a GC, such wide binaries have large 
cross-sections for stellar encounters which can disrupt them, eject them from the cluster core, or 
after multiple collisions, induce significant eccentricities (36) in their initially circular orbits (31}. 
For these pulsars, the measured eccentricities are e~0.02, significantly larger than predicted for 
binary MSPs with helium WD companions in the Galaxy (31 ). If they are located near the cluster 
center, where the stellar densities are >10^pc~^ (E5l) . the time scale for interactions to produce 
such eccentricities is 10*^—10^ yr S36[ . consistent with the >10^ yr ages of most cluster pulsars. 
However, the 60-day orbit of Ter 5 E implies a time scale near the low end of that range, possibly 
too short for MSP lifetimes. This may indicate that it resides further from the center of Terzan 5: 
> 10^ yr spent in the core could result in many interactions which would either destroy the binary 
or induce much larger eccentricities. The eccentricities of Ter 5 I and J are too large to have been 
produced by this method, especially given the relatively compact nature of their orbits, indicating 
that they were formed by a different mechanism. 

At least two of the new binary pulsars, Ter 5 O and P, have been observed to eclipse, al- 
though given the factor of 10 difference in the inferred companion masses (m2^O.O36M0 for 
Ter 5 O and >0.38 M© for Ter 5 P) and eclipse durations (~0.05Porb for Ter 5 O and >0.5Porb 
for Ter 5 P), the systems are unalike. Ter 5 O, the third fastest MSP known, is similar to the 
Black- Widow eclipsing MSP B 1957-1-20 (1571) with a very low-mass companion and an eclipse 
duration (~15— 20min) corresponding to a physical size (~0.03 Rq) much larger than the com- 
panion star's Roche lobe. Such systems are common in GCs (0^ . Ter 5 P, the fourth fastest 
MSP known, is unusual and most similar to the binary MSP NGC 6397 A, whose peculiar red 
straggler or sub-subgiant companion causes irregular and long-duration eclipses of the radio 
pulses @. Ter 5 P's eclipses are also irregular and correspond to an eclipsing region of several 
solar radii in size, indicating the companion is likely a peculiar evolved star as well. If so, the 
system may have been created after at least one exchange encounter where the MSP's original 
companion (which had spun up the pulsar) was ejected and replaced by a main sequence star. 
Such an encounter might eject the system from the core of the cluster as has been observed for 
NGC 6397 A [although alternative formation scenarios for this pulsar have been proposed as 
well (59)]. 

The pulsars Ter 5 I and J are in highly eccentric orbits with companions of at least 0.24 and 
0.38 M0, respectively (Fig. 2). The companions are not sufficiently massive to be NSs unless the 
systems have improbably low inclination angles i (i.e. are almost face-on). If they were main- 
sequence or giant stars, eclipses would likely be observed for a variety of inclinations given the 
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compactness of the orbits (orbital separations of ~5Ro). Additionally, the orbital circularization 
time scales would be of order 10^ yr for both pulsars and so no eccentricity should currently 
be observed I^Q. The companions are therefore probably WDs. However, the pulsar recycling 
scenario generates MSP— WD binaries in nearly circular orbits due to tidal interactions during 
mass-transfer and therefore did not produce these eccentric systems. A possible scenario involves 
the off-center collision of a NS with a red giant which would disrupt the giant's envelope and 
leave the core (now the WD) in an eccentric orbit about the NS (4l\i. An alternative scenario 
involves the exchange of an initially isolated WD into a binary consisting of an MSP and the low- 
mass WD that recycled it. However, the observed pulsars are only mildly recycled. Furthermore, 
the multiple exchanges required in this scenario imply that the system was created near the core, 
where due to mass segregation the isolated WDs would be more massive than the observed 
companions to Ter 5 I and J. We therefore regard the collision scenario as more likely. 

The high eccentricities of both orbits permitted us to measure their advances of the angle of 
periastron, to = 0.255±0.001° yr^^ for I and 0.327 ±0.004° yr^^ for J, which are likely due to 
general-relativistic modification of the Keplerian elliptical orbit. Under this assumption, the total 
mass can be derived for each system: 2. 17±0.02 Mq for I and 2.20±0.04 Mq for J. Together with 
the measured Keplerian mass functions, in both cases this implies that the most likely NS masses 
are >1.7M0 (Fig. 2). However, classical contributions to lu from tidal or rotational quadrupole 
mass moments in the companion star must also be considered (|?2ll?3b . Tidal deformation is 
insignificant for a WD companion (1?^ . but rotationally induced quadrupoles are possible if the 
WD is rapidly rotating, as there is no reason to expect the spin axis of the companion to be 
aligned with the orbital angular momentum. 

The predicted contributions to cu for Ter 5 I and J due to a rotationally induced quadrupole are 
i^rot~0.01— 0.02° yr^^ times angular and stellar structure factors (\42\[43\l44^ . One of the stellar 
structure factors, ag, can range up to 15 for WDs {42), making ujj-ot a potentially significant 
contribution to the measured lu. A significant tUrot, however, will produce changes in i and hence 
in the projected semi-major axis x = aism{i)/c. The size of this effect may be written as 
^vot ~ ivot/x times a trigonometric factor usually of order unity that depends on i, the angle 
between the WD rotation axis and the orbital angular momentum vector, and the phase of the 
orbital precession For random choices of these angles the trigonometric factor is <10, 

~80% of the time. By incorporating detections from Parkes search observations taken in 1998 
and 2000 {21 ) using the now known orbital ephemerides, we have upper limits on x which imply 
that cjrot is ^0.003° yr^^ times the trigonometric factor for both systems. The magnitudes of these 
rotational-quadrupole contributions are comparable to our current measurement uncertainties for 
CO. We emphasize that this result does not depend on the type of companion and is equally 
valid for main-sequence stars and WDs. Therefore, unless the orientations of the orbital and 
rotational angular momentum vectors are fine-tuned for both systems, the cu measurements are 



6 



well-described by general relativity. Under this assumption, one or both of these NSs are likely 
to be unprecedentedly massive: calculation of the joint probabilities (Fig. 2) indicates that at 
least one of the pulsars is more massive than 1.48, 1.68, or 1.74 Mq at 99%, 95% and 90% 
confidence levels, respectively. A strict upper limit to the masses of both pulsars is 1.96 M©. 
Similar but slightly less stringent evidence from pulsar timing for a massive NS in the Galactic 
low-mass WD— MSP binary J0751 + 1807 has been presented {45). It is an intriguing question 
whether the large masses for Ter 5 I and J resulted from the collision formation mechanism, 
possibly during a short period of hypercritical accretion which partially or completely spun-up 
the pulsars (l?7T) . GBT observations over the next 1—3 yr will measure the varying delays due to 
gravitational redshift and time dilation as Ter 5 I moves in its orbit, which together are known as 
7 (the measurement of 7 for Ter 5 J will take 5 — lOyr due to the relatively slow spin period of 
the pulsar). This parameter, along with general relativity and the already well measured cu, will 
provide a precise mass for the pulsar and will likely rule out several soft equations of state for 
matter at nuclear densities (i4<5i) . 
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Table 1 : Known pulsars in Terzan 5 (17^ . Pulsars listed without orbital parameters are likely 
isolated systems while those marked with an are eclipsing systems. The errors on the dispersion 
measures (DMs) range from 0.01—0.1 pccm^'^ and the errors on the measured flux densities are 
~30%. The flux densities for the eclipsing pulsars include only the times when the pulsar is 
not eclipsed. The light travel time across the projected pulsar semi-major axis is defined as 
X = tti sin(z) / c. Eccentricities listed as "0" are too small to measure at present and have been set 
to zero for orbital parameter fitting. The minimum companion mass m2 was calculated assuming 
a pulsar mass mi of 1 .4 Mq and i=90° except for Ter 5 I and J (Fig. 2). All measured parameters 
were determined using the TEMPO software package {47j. 
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Figure 1: The 1950 MHz GBT+Spigot pulse profiles for each of the pulsars known in Terzan 5. 
All but Ter 5 A, C, and D are newly discovered dT^ . Each profile is the weighted average of 
the best detections of that pulsar and is a measure of the relative flux density as a function of 
rotational phase. Asterisks indicate that the pulsar is the member of a binary system, and the 
length of the horizontal error bar (0.3 ms) is the effective system time resolution. 
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Figure 2: Pulsar mass vs. companion mass diagrams for the two highly eccentric binary pulsars 
Ter 5 I (left) and Ter 5 J (right). The hatched regions are excluded due to the definition of the 
Keplerian mass function (i.e. sin i < 1). The diagonal band in the center of the figure shows the 
total system mass with l-a confidence intervals as measured by the general-relativistic advance 
of periastron co'=0.255±0.001° yr^^ for I and 0.327 ±0.004° yr^^ for J. The marginal distributions 
show the masses of the pulsars and companions assuming a random distribution of inclinations 
(i.e., with probability density flat in cosi). The solid curves in the main plot indicate inclinations 
of (from top to bottom) 1 1.4°, 18.2°, 30°, 60°, and 90°. The grey vertical band shows the range 
of precisely measured NS masses from relativistic binary radio pulsars {46). In both cases the 
median pulsar mass (indicated by the vertical line in the marginal distribution) lies significantly 
above 1.7 M©, implying that one or both of these pulsars is considerably more massive than the 
NSs that have been well measured to date. A strict upper limit to the masses of both pulsars is 
1.96 M«. 
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